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Abstract 

We present the observations of the reflection nebulae IC4954 and IC4955 region 
with the Infrared Camera (IRC) and the Far-Infrared Surveyor (FIS) on board the 
infrared astronomical satellite AKARI during its performance verification phase. We 
obtained 7 band images from 7 to 160 fim with higher spatial resolution and higher 
sensitivities than previous observations. Among the 7 band images, the Sll (11 /im) 
data provide the highest spatial-resolution and most sensitive image, in which 10 
point-like sources are newly detected in the mid-infrared and IC4955 is resolved into 
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two condensations. The mid-infrared images show distinct arc-hke structures, inside 
of which bright sources in the L18W image are located. The mid-infrared color 
of the S9W (9 /im) and L18W (18 /im) bands shows a systematic variation around 
the exciting sources. It is red in the vicinity of the exciting sources, whereas it 
becomes blue in the surrounding regions. This can be interpreted in terms of the 
decreasing contribution of thermal emission from dust grains in equilibrium with the 
incident radiation field as the distance from the exciting source increases. The spatial 
variation in the mid-infrared color suggests that the star- formation in IC4954/4955 
is progressing from south-west to north-east. The FIS data also clearly resolve two 
nebulae for the first time in the far-infrared. The FIS 4-band data from 65 fim to 
160 /xm allow us to correctly estimate the total infrared luminosity from the region, 
which is about one sixth of the energy emitted from the existing stellar sources. There 
is little possibility for the presence of embedded massive stars that have escaped 
detection. Five candidates for young stellar objects have been detected as point 
sources for the first time in the 11 fim image. They are located in the red S9W to 
L18W color regions, suggesting that current star-formation has been triggered by 
previous star-formation activities. A wide area map of the size of about 1° x 1° 
around the IC4954/4955 region was created from the AKARI mid-infrared all-sky 
survey data. Together with the Hi 21cm data, it suggests a large hollow structure of 
a degree scale, on whose edge the IC4954/4955 region has been created, indicating 
star formation over three generations in largely different spatial scales. 

Key words: infrared: ISM — open cluster: individual (Roslund 4) — nebulae: 
individual (IC4954 and IC4955) — stars: star formation 

1. Introduction 

The reflection nebulae IC4954 and IC4955 are located in the Vulpecula constellation on 
the Galactic plane around (/,6) = (66.96°, — 1.26°). These nebulae are associated with the young 
open cluster Roslund 4 (Roslund 1960). The heliocentric distance and the age of this cluster 
are estimated based on the analysis of isochrones as lOMyr and 2.9 kpc by Racine (1996), 
4Myr and 2 kpc by Phelps (2003) (hereafter P03), and 15Myr and 2 kpc by Delgado et al. 
(2004) (hereafter D04). These studies suggest the presence of stars of a relatively wide age 
range and on-going star-formation in the IC4954/4955 region. In this paper we assume the 
distance to be 2 kpc. P03 and D04 also detected three active regions in [SII] images and some 
of them are classifled as Herbig Haro objects. The Ha, [Nil] and [SII] emission lines from the 
cluster members are detected in optical spectroscopic observations (D04). This region is an 
interesting place for the study of star-formation process because of the co-existence of both 
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young and relatively old populations (P03; D04). Except for the dedicated observation of D04, 
the parallax, radial velocity, proper motions and spectral type of the cluster members (or clues 
for the identification of the members) are largely unavailable in existing catalogs because of its 
far heliocentric distance. 

^^CO observations show that at least one cloud of 10^ Mq is associated with this region 
(Leisawitz et al. 1989). In the infrared, these nebulae are observed in the IRAS and MSX 
Galactic surveys. The IRAS Point Source Catalog (PSC) has entries of three point sources in 
this region, whereas the MSX map shows distinguished structures in four mid-infrared bands 
of A(8.28/im), C(11.2yum), D(14.3/im) and E(21.34/im). However, the poor spatial resolution 
of IRAS and the insufficient sensitivity of the MSX 21 fim prevent us from making detailed 
analysis of this interesting region in the infrared. This region is not included in the Galactic 
Legacy Infrared Mid-Plane Survey Extraordinaire with Spitzer (GLIMPSE; Benjamin et al. 
2003). 

AKARI is the first Japanese infrared astronomical satellite dedicated for infrared as- 
tronomy (Murakami et al. 2007). It has two scientific instruments: the Infrared Camera (IRC; 
Onaka et al. 2007b) that covers near- and mid-infrared wavelengths of 2-26 /xm, and the Far- 
Infrared Surveyor (FIS; Kawada et al. 2007) that covers the far-infrared spectral range of 
50-200 /im. The performance of both instruments has been confirmed by observations during 
the performance verification (PV) phase from 2006 April 24 to May 8. The present paper 
reports the results of observations of the IC4954/4955 region with both IRC and FIS taken 
mainly during the PV phase. 

The observations and data reduction are described in §2. The observational results are 
shown in §3. The nature and the origin of the Roslund 4 region are discussed in § 4. Finally, 
we summarize the results in §5. 

2. Observations and Data reductions 

AKARI seven band images were taken by two scientific instruments with three kinds of 
operation mode. The observational parameters are summarized in Table 1. 

2.1. IRC All-Sky Survey (9 and 18 fim) 

The S9W (9 /im) and L18W (18 /im) data of the IC4954/4955 region were taken as part 
of the all-sky survey observations. Both data were taken simultaneously with different channels 
of the IRC, MIR-S and MIR-L, which observe sky positions separated by about 20' in the cross- 
scan direction (Onaka et al. 2007b). In the all-sky survey, the IRC is operated in the scan mode 
(Ishihara et al. 2006a), in which the data of only two rows in the detector arrays are taken with 
the cross-scan width of about 10'. The scan speed is about 215" s~^. The output signals of 
every four adjacent pixels are binned together to reduce the data down- link rate and the virtual 
pixel scale is 9. "36 x 9. "36 in the survey mode. The seconds confirmation is enabled by the 
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Table 1. Observational parameters for IC4954/4955 by AKARI 
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independent data sets taken by the two separated rows, which allow to efficiently reject high- 
energy ionization particle events and largely improves the reliability of the source detection. 
Finer spatial information as well as improved sensitivity can be expected in the data reduction 
process because two rows are sampled and binned in an interlaced manner and the observed 
strip is shifted by 4' in adjacent scan paths with a given region being observed more than twice 
on average. The reset interval is set as 13.5 s, which corresponds to 306 samplings or 48. '9 in 
the scanning direction. The length of a reset in the all-sky survey is set as 2.2 ms (hereafter 
long reset). The IC4954/4955 region was observed during 2006 May 1-8 in the PV phase with 
the descending path of dp/dt<0 and 2006 November 1-8 with the ascending path of df3 /dt> 0, 
where (3 is the ecliptic latitude. Both data are averaged with median in the present analysis. 
The difference between the two data sets is less than 10% in the sky brightness. 

Fragments of the data sets covering IC4954/4955 are retrieved from the IRC all-sky 
survey database and reduced automatically by the pipeline software to create a map (Ishihara 
et al. 2007). The pipeline process includes linearity correction, flat correction, reset anomaly 
correction (see below), source extraction, rejection of high-energy particle events by the seconds 
con&mation, position determination using identified objects, coaddition of the images obtained 
by the two rows, and map construction. 

The flux calibration of point sources in the all-sky survey data are carried out based 
on a large number of detections of hundreds of stars in the standard star networks (Cohen et 
al. 1999; Cohen et al. 2003; Ishihara et al. 2006b). The flux accuracy is estimated as ~ 7% 
for the S9W band and ~ 15% for the L18W band for point sources at the present calibration 
stage. The detection limits (50") for point sources are estimated as 50mJy for the S9W band 
and 120 mJy for the L18W band. The calibration for diffuse emission is on-going and a relative 
error of 30% is assigned at present. 

Positions of the detected objects were ffist estimated based on the data of the gyroscopes 
and the star trackers on board and then refined by using the positions of bright (K<8 mag) 
2MASS PSC sources. The position accuracy is estimated to be 5" at present. 
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2.2. IRC Slow Scan (llfim) 

The Sll (11 /iin) band image was obtained during the PV phase with the IRC slow scan 
mode of the Astronomical Observation Template (AOT) IRC51. IRC51 is a set of round-trip 
scans designed for mapping of a small area of up to ~ 10' x 120' in about 15min of a pointed 
observation (Murakami et al. 2007; Onaka et al. 2007b). A higher sensitivity than in the all- 
sky survey and a wider sky coverage than in imaging observations can be achieved in this 
observation mode. It also allows to observe very bright sources, such as IC4954/4955, which 
will be saturated in the imaging observation. The scan speed was set as 30"s~^ and the Sll 
filter was selected. The observation of IC4954/4955 was carried out on 2006 May 1. 

The focal plane arrays were operated in the same scan mode as in the all-sky survey 
observations except that the data binning is not made and the full spatial information is ob- 
tained in the cross-scan direction. Confirmation of source detection was made by round-trip 
scans. The grid size of slow scan observations is set as l."32 x 2. "34. It is determined by the 
physical pixel size in the cross-scan direction (2. "34), and by the sampling rate and the scan 
speed in the scanning direction (l."32). It oversamples the image size of the imaging mode 
at 11 /im (~ 4. "8; Onaka et al. 2007b). The time spent for a reset was shortened to 68 /is 
(short reset). The charge integration curve shows an anomalous behavior immediately after 
a reset (reset anomaly). The magnitude of this effect depends on the time spent for a reset. 
We confirm that the short reset significantly improves the reset anomaly effect. Gaps in the 
observed area due to the reset was thus reduced thanks to the combination of the slower scan 
speed and the short reset. The reset period was also shortened to 2.24 s because the gap due 
to the reset could be ignored. Consequently the effect of saturation was also reduced. There is 
a drift in the offset level owing to the temperature drift (Ishihara et al. 2003). It is corrected 
by referring to the signal level during the reset. The signal during the reset corresponds to 
the output with the input blocked, and thus indicates the offset level of the preamplifier. The 
detection limit is determined by the read-out noise of the detector rather than the zodiacal 
background fluctuation because of the short sampling time. 

The data reduction largely follows that for the all-sky survey except for a few parameters 
adjusted for the finer pixels, such as fiat correction, and a custom-designed map reconstruction 
module is developed. At first the data taken in a single scan are arranged into a 10' x 60' strip 
image by assuming that the sampling rate and the scan speed of the satellite are constant. 
Individual sets of images obtained by two separated rows are aligned with each other and 
coadded into a combined image. The position reconstruction of the map is made based on the 
association with detected 2MASS sources as in the all-sky survey data. The positional accuracy 
of the data is estimated to be 5" at present. 

The in-fiight calibration of the IRC has been carried out for all the imaging filters in 
pointed imaging observations based on observations of the standard star networks (Onaka et 
al. 2007b). For the all-sky observations and slow-scan observations, however, most observa- 
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tions have been made only with the S9W and L18W bands. The flux calibration for the Sll 
(11 /im) band image in the slow-scan observations is performed in an indirect way because no 
observations of standard stars have been made in the same scan speed with the same filter. 
First, the calibration for the S9W in the slow-scan mode with the scan speed of 30" s~^ is made 
by the comparison of fluxes derived from slow-scan observations with those from the all-sky 
survey for the same stars. Then the calibration for the Sll band with the scan speed of 30 " s~^ 
is estimated by assuming the relative calibration between the S9W and Sll bands in pointed 
imaging observations. Taking account of the uncertainties in the indirect calibration, we set a 
conservative uncertainty of 30% in the Sll data at present. 

2.3. FIS Slow Scan (65,90,140,160 fim band) 

FIS observations of IC4954/4955 were executed with N60 (65 /xm), WIDE-S (90 /im), 
WIDE-L (140 /zm), and N160 (160 /im) bands in the 2-round-trip slow-scan mode (FISOl) with 
a scan speed of 15"s~^ on 2006 May 3. The four-band data were taken simultaneously with 
correlated double sampling (CDS) in about 15min of a pointed observation. Details of the 
observation scheme and the data reduction of the FIS slow scan data are described in Kaneda 
et al. (2007b) and Suzuki et al. (2007). At present, the systematic flux calibration errors for 
the CDS mode are estimated to be 30, 40, 50, and 50% for N60, WIDE-S (90 /im), WIDE-L 
(140 /xm), and N160, respectively. The relative position accuracy among the FIS bands is better 
than 10" because they were taken simultaneously, but the absolute accuracy is estimated to 
be about 1'. Thus all the FIS images are aligned relative to the IRC images in the equatorial 
coordinates. Since the IRC data were taken at a different epoch from that of FIS, the relative 
position accuracy is estimated to be about 30" at present. 

3. Results 

3.1. Infrared images of I C4 9 54 and I C 49 5 5 

The images of IC4954/4955 in AKARI S9W, Sll, L18W, N60, WIDE-S (90 /xm), WIDE- 
L (140 /im), and N160 bands are shown in figure 1. The northern nebula in the image is IC4955 
(enclosed by the green line in figure la) and the southern nebula is IC4954 (enclosed by the 
red line). The AKARI images reveal several distinct characteristics of the nebulae. 

The AKARI S9W map by the all-sky survey is created with a grid size of l."56 x l."56 
from the virtual pixel size of 9. "36 x 9. "36. It is in good agreement with the MSX A-band 
(8.28 /im) image. The AKARI S9W band has a better spatial resolution and the relative spectral 
response similar to that of MSX A-band except for the inclusion of the 11.2 /xm unidentified 
infrared (UIR) band (Onaka et al. 2007b). Three IRAS point sources exist in this region as 
indicated in the S9W image. The location of IRAS 20026+2906 does not match with the bright 
spot in the S9W images (source F in figure lb) probably due to the complicated background 
structures in this region. 
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The Sll image has the finest resolution and the deepest sensitivity among the seven 
bands owing to its observing mode. Point-hke objects are clearly resolved and denoted by A-L 
(for details see §3.2). Distinct extended structures are denoted hj a — h. The northern nebula 
IC4955 is spatially resolved into two condensations (a and d) and two arc-like structures {b and 
c). The two arcs may overlap with each other on the line-of-sight . Arc-like structures are also 
found in IC4954 (e and /) . All the arcs point toward the north-east of the nebulae and extend 
over the outer edge of the nebula emission in the optical image. 

The AKARI L18W image has a higher signal-to-noise ratio and higher spatial resolution 
than the MSX E-band (21.2 yum) image. Point-like objects (g-n) are clearly seen in figure Ic. 
The spatial distribution of the 18 fim emission seems to be significantly different to that of the 
S9W or Sll emission. See §4.1 for the discussion on the morphologies of the S9W and L18W 
images. 

The AKARI N60 (65 /im) and WIDE-S (90 /im) image surpass the IRAS 60 or 100 /im 
images. The two nebulae are spatially resolved in this wavelength range for the first time. The 
faint arc-like structure toward the south-east is also detected. The WIDE-L (140 /im) and N160 
(160 /im) data are unique to AKARI. These two bands make significant constraints on the dust 
emission (figure 2) for the estimate of the total infrared luminosity (§ 4.1). 

The total fiuxes from IC4954 and IC4955 are derived from the integrated signals in the 
regions enclosed by the red and green lines in figure la and the sky background is estimated 
from the surrounding region and subtracted. The results of IC4954 and IC4955 are summarized 
in Table 2 and their spectral energy distributions (SEDs) are plotted in figure 2. The uncer- 
tainties in Table 2 include the systematic errors. The MSX image data at band A (8.28 /im) 
C (12.13/im), D (14.65 /im), and E (21.3/im) are also integrated over the same regions and 
the sky background is subtracted similarly. The results are plotted in figure 2. Except for 
the band D data, the IRC data are in fair agreement with the MSX data, taking account of 
the uncertainties and the differences in the band profiles. The band D data are fainter even 
compared to other MSX band data for both nebulae. 

In figure 2 also plotted are the IRAS PSC data. It is not possible to estimate the fiuxes 
for IC4954 and IC4955 directly from the IRAS data since the sources are not clearly resolved. 
In the figure, the fiuxes of IRAS 20028+2903 are plotted as IC4954, whereas the sum of the 
fiuxes of IRAS0026+2906 and IRAS2007+2905 are indicated as IC4955. Except for the 60 /im 
data, the IRAS fiuxes are lower than the AKARI/IRC and FIS fiuxes and the agreement is not 
very good. The difference can be attributed to the fact that the IRAS data do not include the 
diffuse emission correctly in addition to the position mismatch of IRAS 20026+2906. 

3.2. Newly detected point-like sources 

The fluxes and positions of 12 point-like objects detected in the Sll image (indicated by 
A-L in flgure lb) are summarized in Table 3. They all have corresponding 2MASS sources as 



7 



Table 2. Photometric results of IC4954 and IC4955* 



AKARI/IRC 
Object 


results 
S9W 
9 jim 


Sll 
11 fini 


L18W 
18 /im 


N60 
65 /im 


WIDE-S 
90 /im 


WIDE-L 
140 /im 


N160 
160 /im 


IC4954 
IC4955 


62.2+18.7 
33.6+10.1 


61.1+18.3 
33.4+10.2 


56.2+16.9 
31.2+9.36 


790+240 
330+100 


860+260 
340+100 


2500+1000 
1700+680 


1890+570 
850+260 


MSX data 
name 


8.28 fj.m 


12.13 /im 


14.65 /im 


21.34 /im 








IC4954 
IC4955 


55.15+0.26 
28.86+1.44 


75.15+2.25 
40.25+1.21 


32.29+1.29 
16.10+0.64 


66.84+4.01 
36.21+2.17 








IRAS PSC sources 
IRAS name 12 /im 


25 fim 


60 /im 


100 /im 








20026+2906 
20027+2905 
20028+2903 


2.741+0.66 
3.356+0.30 
20.26+0.81 


<0.25 
14.06+0.70 
47.29+1.89 


<3.186 
<540.3 
540.3+48.6 


<66.76 
<50.71 
1177+129 









* Units are in Jy. 



indicated in the table. The SED of each source is plotted in figure 3, including the 2MASS data. 
There is a MSX source (G066. 9971-01. 2247) close to source F, whose C band flux (12.13 /im) is 
much larger than the Sll flux (see Table 3), suggesting that the MSX flux includes contributions 
from the surrounding diffuse emission since the source is located in an extended emission region. 
In the vicinity of source K, there are two MSX sources (G066. 9646-01. 2783 and G066.9609- 
01.2776), whose C band fluxes are also very large (see Table 3). G066. 9646-01. 2783 seems to 
correspond to source K in position and G066. 9609-01. 2776 seems to be located close to another 
optically bright source, 2MASS 20045331+2911469, around which the Sll image does not detect 
any point source. The large MSX fluxes can also be attributed to the diffuse emission around 
the source. Except for F and K, 10 sources are detected as point sources for the flrst time in 
the mid-infrared by the present IRC observations. 

4. Discussion 

4.I. Infrared characteristics of IC4954 and IC4955 

The total infrared luminosity from the nebular region can be estimated from the AKARI 
observations. The AKARI IRC and FIS data (figure 2) are fitted by a three-temperature dust 
model with the emissivity proportional to A~^. The three temperatures are fixed as 20, 40, and 
230 K for both nebulae. The fitted models are shown by the solid and dotted lines in figure 2 
for IC4954 and IC4955, respectively. The total infrared luminosity Lin is then calculated by 
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Table 3. Point sources detected in the Sll image 



ID 


R.A. Decl. 


2MASS ID 


1 1 fiia flux 


MSX source 


MSX band C flux 




(J2000.0) 




(mJy) 




(Jy) 



A 


20 


04 


56 


+29 


11 


17 


2004568H 


h2911154 


6.7+0.3 






B 


20 


04 


46 


+29 


11 


16 


20044653H 


h2911167 


19.6+0.3 






C 


20 


04 


49 


+29 


11 


48 


20044974H 


h2911486 


55.6+1.2 






D 


20 


04 


43 


+29 


12 


49 


20044366H 


h2912496 


13.1+0.3 






E 


20 


04 


45 


+29 


13 


39 


20044585H 


h2913389 


22.3+0.5 






F 


20 


04 


45 


+29 


15 


14 


20044568H 


h2915152 


23.0+0.8 


G066.9971-01.2247 


3.4441 


G 


20 


04 


44 


+29 


13 


50 


20044473H 


h2913512 


7.7+0.2 






H 


20 


04 


46 


+29 


10 


21 


20044676H 


h2910221 


5.2+0.2 






I 


20 


04 


47 


+29 


09 


59 


20044710H 


h2910006 


39.2+0.4 






J 


20 


04 


50 


+29 


11 


05 


20045044H 


h2911060 


45.7+1.1 






K 


20 


04 


52 


+29 


11 


43 


20045278H 


h2911435 


208.5+4.4 


G066.9646-01.2783* 


5.4393 






















G066.9609-01.2776t 


3.6187 


L 


20 


04 


54 


+29 


12 


07 


20045478+2912080 


14.3+0.3 







* Most probable MSX source that corresponds to K. 

^ No corresponding point source in the Sll image (see text). 



1=1'' 



where Tj = 20,40, and 230 K for i = 1,2, and 3, respectively, and Ci is the fitting constants. We 
obtain Lj^ = 9.5 x 10^ and 4.2 x IO^Lq for IC4954 and IC4955, respectively, the sum of the two 
nebulae being 1.4 x IO^Lq. 

The stellar luminosity available for dust heating can be estimated by summing up the 
stellar luminosity of all the stars located in the region as 

-^stellar = ^ -^*- 
members 

The luminosity of each star is estimated simply based on its spectral type with the 
assumption that they are on the main sequence (Schmidt-Kaler 1982; Johnson 1996; De Jager 
et al. 1987). We select 13 stars that have known spectral types and are assigned to the member 
of this region (D04). We also select other 11 stars whose spectral types are estimated to be 
earlier than B7 based on their color. A total of 24 stars in this region are included in the 
estimate of /^stellar- The spectral type is taken from the spectroscopic classification or estimated 
from the reddening corrected U — B and B — V m D04. The membership assignment and the 
UBVRI photometric results (D04) were obtained from WEBDA^. The location of the 24 stars 
used in the estimation is shown by the crosses in figure 4a. The 08en star reported in Merrill 
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et al. (1942) located in the south-west of the nebulae (shown in figure 4b) is not included in 
the estimation because it seems to be too far from the nebula region. Even if we include this 
star, the total stellar luminosity will be increased only by 15%, and it will not affect the present 
conclusion. The total stellar luminosity in this region thus estimated is 8.3 x W^Lq. It is about 
6 times of Ljr. Leisawitz & Hauser (1988) have indicated that typically only a small fraction 
(~ 0.2) of the cluster luminosity is absorbed by dust grains. The present result is in good 
agreement with their results, suggesting that it is not very likely that there are very luminous 
stars hidden in dusty environments in this region. 

With the dust mass emissivity at 100 /xm being 0.6gcm~^ (Hildebrand 1983), the dust 
mass associated with IC4954 and IC4955 is derived to be about 40 and 20 Mq, respectively. 
The typical size of the nebulae is estimated from the intensity contours of 10% of the peak to be 
about 1.'3 or 2 X 10^® cm. Assuming that the gas-to-dust ratio is 100, the average gas density of 
the infrared emitting medium is estimated to be about (2 — 4) x 10^ cm~^. This is in the range 
of the density of dense photo-dissociation regions (PDRs), such as the Orion region (Tielens & 
HoUenbach 1985), indicating that the IC4954/4955 region is still a young star- forming region 
rather than diffuse PDRs, such as the Carina nebula (Mizutani et al. 2004). 

The color images of IC4954/4955 in the optical, mid- infrared, and far- infrared are shown 
in figure 4. In the mid-infrared, the spatial distribution of the 9 /xm emission is similar to that 
of the 11 /im emission, whereas the difference in the spatial distribution between S9W and 
L18W images is remarkable. Figure 5 shows the relative spectral response of the S9W, Sll, 
and L18W bands. Also plotted by the thin solid line is an interstellar cirrus spectrum taken 
with Spitzer/IKS for a reference. The S9W band includes the major UIR bands at 6.2, 7.7, 8.6, 
and 11.2 /im except for the 12.7 fim band, and is not affected by the continuum emission longer 
than 12 /im. The Sll band is affected by only part of the UIR 7.7/im, but includes the UIR 
11.2 and 12.7 /im, and the [Nell] 12.8 /im line emission. In the L18W the UIR 17 fim complex is 
included and the main contributor is the continuum emission longer than 15 /im. No appreciable 
difference seen between the S9W and Sll images suggests that the spectrum shape between 
6 to 13 /im does not vary significantly in the region and that the line emissions from ionized 
species is insignificant compared to the UIR band emission. The difference between S9W and 
L18W should be related to the variation in the continuum emission longer than 15 fim. The 
S9W image shows arc-like structures clearly and the L18W image indicates high concentration 
of the emission in narrow areas. The difference can be seen noticeably in the mid-infrared color 
image (figure 4b). The emission in L18W is stronger inside the arcs as indicated by red. Most 
of the regions highlighted by red color are associated with B-type stars, which are supposed 
to play as heating sources of the region: in the region g, there are a B and a B9 star; the 
region h is associated with a B2 and a B7 star; a B-type star in the region n is suggested to 
not belong to the member of this region on the basis of the radial velocity (D04) and thus 
may not be related to this region. The other nearby B5 star must be a heating source of this 
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region. The red regions match also with bright regions in the Ha image (D04). We surmise that 
the red regions are directly heated by B-type stars and probably associated with ionized gas, 
whereas the arcs represented by white color in figure 4b are PDRs surrounding them, which 
are characterized by the strong UIR emissions in the S9W (cf., Onaka 2004). The strength 
of the incident radiation field at the arcs {b, c, and /) is estimated based on the projected 
distance from these heating sources as 900, 3000, and 800, respectively, in units of the solar 
neighborhood value (1.6 x 10~^Wm~^, Habing 1968). These are comparable with those of 
typical PDRs (cf., Mizutani et al. 2004). 

IRAS observations indicate the increase of the ratio of 25 fim to 12 /im intensities in 
the vicinity of heating sources for several cases (Boulanger et al. 1988; Ryter et al. 1987). 
The UIR band emission largely contributes to the IRAS 12 fim band. The variation in the 
intensity ratio has been attributed to the increasing contribution from thermal emission of 
large grains rather than the destruction of the UIR band carriers. Recent investigations on the 
infrared diffuse radiation of our Galaxy and the Large Magellanic Cloud have in fact shown 
that the variations in the infrared SED of the diffuse emission can reasonably be interpreted in 
terms of a superposition of the model emissions, in which the contribution from dust grains in 
radiative equilibrium becomes large in the vicinity of heating sources together with the effect 
of destruction of the UIR band carriers in strong radiation fields (Sakon et al. 2006; Onaka 
et al. 2007a). The band profile of the IRC S9W is shifted to shorter wavelengths compared 
to the IRAS 12 fim band and thus S9W probes the UIR band emission more sensitively and 
is less subject to thermal emission in longer wavelengths (Onaka et al. 2007b). Consequently 
the effect of thermal emission near the exciting source appears clearly in the color map of 
S9W/S11/L18W. The red color in figure 4 points to the region strongly heated by heating 
sources, whereas the white color indicates the region where the UIR band emission is dominant. 
The observed color variation is well accounted for by the increasing contribution from thermal 
dust emission. It further indicates that the UIR-band dominating infrared bright regions (white 
in figure 4b) are always located in the north-east side, suggesting the presence of high density 
regions in this side. 

To examine the distribution of young stellar object (YSO) candidates in this region, ob- 
jects with red color in the near-infrared are selected from the 2MASS catalog. The interstellar 
reddening to the IC4954/4955 region is estimated to be E{B — V)= 0.91 (D04), which corre- 
sponds to [H — K) = 0.18 (Rieke & Lebofsky 1985). Taking account of the internal extinction, 
we conservatively set the criteria that H — K > 0.9 and J — H > 0.9 to select YSO candidates. 
They are shown by the red circles and other blue 2MASS sources are indicated by the blue 
circles in figure 4b. Most of them are located in the white color region as well as in between 
the two nebulae, where the optical image shows few stars (see below). Only a few red objects 
are present inside the arcs, indicating that YSO candidates are concentrated in dense regions 
slightly away from those directly heated by B-type stars. Five of the 2MASS red objects are 
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detected in the Sll image (C, F, J, K, and L). All of them show YSO-like SEDs in the near- 
to mid-infrared (§ 3.2; cf., Whitney et al. 2003; Reach et al. 2004), supporting the validity of 
the adopted criteria. There is also one red object in the relatively red color region without 
any corresponding early-type stars. It may also play a role of an embedded heating source in 
the region. A summary of the relations among the various extended [a-n) point-like sources is 
given in Table 4. 





Table 4. 


Structures and point sources in the mid- infrared of the IC4954/4955 region 


Sll 




S18W 


Comments 


Extended Point-like 


Extended 




a 


F 


9 


The peak positions of a and g do not coincide. 








F is located at the peak of a. 


b 




i 


The peak positions do not coincide. 








There may be a hidden source at i. 


c 


E, G 


h 


The peak positions of c and h do not coincide. 








E and G are heating sources for h. 


d 




3 


The peak positions are matched within the uncertainties. 




K 


k,l 


K coincides with k. No counterpart is seen for /. 




A 


n 


The position of A coincides with n. A is listed as a non-member. 




J 


m 


J is a heating source for m. 



Comparison of figures 4a and c indicates that the density of stars is low in between IC4954 
and IC4955 in the optical image, whereas the dust emission is seen in the corresponding region 
in the WID15-L (140 /im) image. The size of the region (120") is sufficiently larger than the 
FWHM of the beam size of the WIDEl-L band (~ 58"; Kawada et al. 2007), thus the presence 
of the far-infrared emission in between the two nebulae is not spurious. It can be surmised that 
the region is a dark nebula filled with dust grains. The region around source F in figure lb also 
appears dark in the optical image, suggesting that source F is surrounded by a large amount 
of dust. 

The mid-infrared color variation, the location of early-type stars, and the distribution of 
red objects strongly suggest that YSOs are being born in the region on the north-east side of 
the arcs. Their formation may be triggered by the existing early-type stars and star-formation 
is propagating from south-west to north-east. This picture is revealed by the high spatial 
resolution (especially in 11 /im) and the multi-band (especially 140 or 160 /xm) infrared data of 
AKARI. Because of the arc-like shape rather than the cometary shape, it is conjectured that 
the regions where YSOs are concentrated in IC4954/4955 are not pre-existing clouds imploded 
and/or evaporated by the effect of stellar winds, such as in the Elephant Trunk Nebula (Reach 
et al. 2004), but are rather triggered by the 'collect and collapse' type mechanism with stellar 
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winds of the heating stars (Elmegreen 1998). Rather uniformly distributed YSOs may also 
support this interpretation (e.g., Efremov & Elmegreen 1998). 

4-2. Origin of Roslund 4 

To investigate the origin of the IC4954/4955, mid-infrared maps of the surrounding 
region of about 1° x 1° at 9 and 18 fim are created from the IRC all-sky survey data. The 
S9W/L18W color image is shown in figure 6a with the contours of the 100 /im data of the 
IRAS Sky Survey Atlas (ISSA, Wheelock et al. 1994). It indicates that there is a cavity of low 
mid- infrared emission centered around {a, 6) =(20h 03m, -|-29°.00) and that the IC4954/4955 
region is located on an edge of the cavity, image (figure 1). 

To examine and confirm the presence of the cavity. Hi 21 cm data of the region are 
obtained from the Canadian Galactic Plane Survey (CGPS)^ and examined. Figure 7 shows a 
plot of the intensity vs. velocity of the Hi data. It indicates that the Hi gas in this region has 
a velocity range of 0-30kms~^. This is in good agreement of the velocity range of CO emission 
associated this region of 6-19 km s~^ (Leisawitz et al. 1989). The Hi intensity integrated over 
0-30kms~^ is shown in figure 6b, which clearly supports the presence of a low density cavity 
around the center of the map. In figure 6b also plotted by the crosses are YSO candidates. They 
are selected from the 2MASS catalog based on the same criteria, H — K > 0.9 and J — H > 0.9, 
as in figure 4b. In addition a condition of 13.5 < K < 15.0 is estimated from the K magnitude 
of the YSO candidates in figure 4b and applied to exclude foreground and background sources. 
The symbols a, /5, and 7 indicate the regions where concentration of YSOs is seen. YSO 
candidates seem to be located on the edge of the cavity. It is most clearly seen in the eastern 
edge, which includes the IC4954/4955 region as well as the regions a and /3, however, only a 
few YSOs are found at the north and south edges. 

The radius of the cavity is estimated as about 10 pc. The cavity is elongated and partly 
collapsed. Such structures are commonly seen in Galactic bubbles (Churchwell et al. 2006). 
The cavity of the IC4954/4955 may have been created by supernovae (SNe) or 0-type stars. 
Progenitor candidates have been searched for in the Tycho-2 spectral type catalog (Wright et 
al. 2003), All-sky compiled catalog (Kharchenko 2001), Mean radial velocities catalog (Barbier- 
Brossat & Figon 2000), and SNRs catalog (Green et al. 2006). but any progenitors of such a 
kind are not found. Any SNR like structures are not found either in the soft X-ray maps of 
ROSAT or in the NVSS 1.4 GHz map (Condon et al. 1998), although there seem to be some 
point-like condensations in the region. The cavity of the size of ~40pc in IC1396 is made 
by a single 06.5 star (Patel et al. 1995). Thus the 08en star located at the south-west of 
the IC4954/4955 (figure 4b) could be the progenitor because the dense region is spread at the 
north-east of the cavity. However the star may be too young if the age of Roslund 4 is 15 Myr 
and the star-formation was delayed relative to the trigger event because of the internal motions 

^ http: / / www.ras.ucalgary.ca/ CGPS 
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(Elmegreen 1998). Since SNRs disappear in the time scale of ~1 Myr and the age of the cavity 
must be older than 4 Myr based on the age of the Roslund 4, it is not unexpected that there is 
no direct sign for the progenitor. The Cepheus bubble (Patel et al. 1998) is also thought to be 
formed by 0-type stars or supernovae, but shows no direct hint of them. The present AKARI 
observations have revealed another bubble of small scale, which seems to trigger past (4-15 Myr 
old) star- format ion in IC4954/4955. It also indicates the 'current star-formation' in the region, 
which is triggered by stars of the second generation. This trigger seems to be different from the 
'globule-squeezing' type in the Elephant Trunk Nebula near the Cepheus bubble. The present 
observations indicate triggered star- format ion over the three generations in a different spatial 
scale in the IC4954/4955. 

5. Summary 

With the two scientific instruments (IRC and FIS) on board AKARI we obtained 7 band 
images of IC4954/4955 from 7 to 160 /xm with higher spatial resolution and higher sensitivities 
than previous observations. Based on these observations we obtained the following results. 

(1) The mid-infrared images at 9, 11, and 18/im reveal several distinct structures in 
the region. Particularly they show three arc-like structures that seem to have been created by 
stellar winds from the existing B-type stars. The difference between the S9W (9 fim) and Sll 
(11 nm) images is not significant, indicating that the variation in the spectrum of 6-13 fim is 
not large and that the contribution of line emission from ionized gas, such as [Nell] 12.8 yum, is 
relatively small compared to the UIR band emissions at 6.2, 7.7, and 11.2 /im. 

(2) The S9W (9 fim) to L18W (18 fim) images appear systematically different from each 
other. The L18W emission is strong near the exciting stars, whereas the S9W emission is 
dominant in the surrounding region. The S9W filter probes the UIR band emission more 
effectively than the IRAS 12 fim because its band profile is shifted to shorter wavelengths. We 
interpret the systematic mid-infrared color variation in terms of the decreasing contribution of 
thermal dust emission with the distance from the exciting source. The color variation clearly 
indicates the location of exciting stars, suggesting that the star-formation in IC4954/4955 is 
progressing from south-west to north-east. 

(3) Young stellar objects are distinguishable for the first time at 11 fim. They are 
located in the large S9W to L18W ratio regions, suggesting that current star- formation has 
been triggered by previous star-formation activities. 

(4) The FIS data allow to correctly estimate the total infrared luminosity from the region, 
which is about one sixth of the energy emitted from the existing early-type stars. Thus there 
is little possibility that embedded luminous stars have escaped detection. It also suggests that 
the total dust mass of the infrared emission is about 60 Mq. The 140 and 160 fim images reveal 
the presence of a high density region between IC4954 and IC4955, which is also supported by 
optical images. 

14 



(5) The IRC all-sky survey data together with the Hi 21 cm data further suggest the 
presence of a bubble- like structure of a degree scale, on whose edge the IC4954/4955 region is 
located, indicating triggered star formation over three generations. 

The IC4954/4955 region is not a massive star-formation region, and is currently pop- 
ulated by B-type stars. The mass of the infrared emitting material is also not huge. The 
suggested hole of the interstellar matter is not large and thus should be created by a less 
energetic event compared to IC1396. The present observations suggest that even in such a 
moderate star- forming region, the sequential star-formation occurs and is on-going at present. 
Medium-scale star-formation could be common and should be studied in future investigations. 

AKARI continues to provide significant data for the study of interstellar medium and 
star forming regions owing to its high sensitivity, wide wavelength coverage and wide mapping 
capability, during the course of its mission. 
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Fig. 1. IC4954/4955 images by AKARI observations: (a) S9W (9/im), (b) SU (11 Mm), (c) L18W 
(18 Aim), (d) N60 (65 ^im), (e) WIDE-S (90 /xm), (f) WIDE-L (140 ^m), and (g) N160 (160 //m) bands 
in the equatorial (J2000) coordinates. The integrated areas in the photometry for each nebula are shown 
on the S9W image (a) with tlic red Unc for IC4954 and the green line for IC4955. Three IRAS point 
sources are also indicated by the plusses. Tlie syiftols a-f in the Sll (b) and g-n in the L18W (c) 
indicate distinct structures in the images. The symbols A-L sliow point-like sources in the Sll image (b). 



10 100 XQQO 

Fig. 2. The 9-160 /zm spectral energy distribution of IC4954 and 1C4955 measured by AKARl The soHd 
bars show the observational data for IC4954 and the dotted bars for IC4955 derived from the AKARI ob- 
servations. The triangles and squares indicate those estimated from the MSX data for 1C4954 and IC4955, 
respectively. The circles and crosses are those from the IRAS PSC data for 1C4954 and IC4955, respectively 
(see text for details). The solid and dotted curves are fitted 3-tempcrature dust models (see also text). 
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Fig. 4. Synthesized color images of IC4954/4955 in the equatorial coordinates: (Left) Optical three- 
-color image of B (blue), R (green) and Z (red) of the DSS image; (Middle) mid-infrared three-color 
image made from the S9W (blue), SIX (green) and L18W (red) images taken with the AKARI/IRC. 
Optically known early-type stars are indicated by the crosses with the spectral types (Delgado et al. 
2004; MerriU et al. 1942). The symbols of C, F, J, K, and L in blue show point-like objects in the Sll 
image, a-f (green) indicate structures seen in the S9W or Sll images, and g-n (yellow) denote struc- 
tures seen in the L18W image as depicted in figures lb and c. Red {H — K > 0.9 and J — H > 0.9) 
sources selected from the 2MASS catalog are shown by the orange circles. The spectral types of the 
stars are indicated in white. (Right) Far-infrared three-color image of IC4954/4955 made from the 
N60 (blue), WIDE-S (90 /im) (green) and WIDE-L (140 yum) (red) images taken with the AKARI/FIS. 
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Fig. 5. Relative band response of the S9W (thick solid line), Sll (thick dashed line), and 
L18W (thick dotted line) in units of electron/energy normalized at the peak. An inter- 
stellar cirrus spectrum at {I, b) ~ (355. °2, 0.°02) taken with Spitzer/IRS is also shown as 
a reference by the thin solid line. It is in units of Jy, but arbitrarily scaled. 
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Fig. 6. (Left) Color image made from the S9W and L18W all-sky survey data of 
AKARI/IRC aromid IC4954/4955. The contours indicate the ISSA 100 ^m inten- 

sity. (Right) CGPS HI 21cm map integrated over the velocity range of 0-30 km s~^ 
(see text). The crosses indicate red sources selected from the 2MASS catalog. 
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Fig. 7. Intensity vs velocity plot of the HI 21cm data (CGPS). The solid line indicates at the 
IC4954/4955 region, whereas the dashed line corresponds to the center of the cavity (figure 6). 
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